The diffuse Galactic 7-rays from dark matter annihilation 
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The diffuse Gaiactic 7-rays from EGRET observation sliows excesses above 1 GeV in comparison 
with the expectations from conventional Galactic cosmic ray (CR) propagation model. In the 
work we try to solve the "GeV excess" problem by dark matter (DM) annihilation in the frame of 
supersymmetry (SUSY). Compared with previous works, there are three aspects improved in this 
work: first, the direction-independent "boost factor" for diffuse 7-rays from dark matter annihilation 
(DMA) is naturally reproduced by taking the DM substructures into account; second, there is no 
need for renormalization of the diffuse 7-ray background produced by CRs; last but not the least, 
in this work our new propagation model can give consistent results of both diffuse 7-rays and 
antiprotons, by directly adding the signals from DMA to the diffuse 7-ray background. This is a 
self-consistent model among several possible scenarios at present, and can be tested or optimized 
by the forthcoming experiments such as GLAST, PAMELA and AMS02. 



The diffuse Galactic 7-rays are produced via interac- 
tion of CRs with the interstellar medium and radiation 
field. However, the spectrum of the diffuse 7 rays mea- 
sured by EGRET shows an excess above I GeV 1] in 
comparison with the prediction based on the conven- 
tional CR model, whose nucleus and electron spectra are 
consistent with the locally observed data. The discrep- 
ancy may indicate large-scale proton or electron spec- 
trum, which determines the diffuse 7-rays, different than 
the local measured one, or the existence of exotic sources 
of diffuse continuum 7-ray emission. 

A harder nucleon spectrum with power-law index of 
—2.4 ~ 2.5 has been proposed in Ref. [2] to solve the 
"GeV excess" problem. However, it has been pointed 
out that such a hard nucleon spectrum will overproduce 
secondary antiprotons and positrons 3] , which has effec- 
tively been excluded by recently high energy p/p ratio 
measurements [1]. A hard electron spectrum is studied 
in Ref. Q while this hypothesis also suffers difficulties, 
e.g. it produced too many 7-rays at higher energies and 
couldn't be compatible with the local electron spectrum 
. For the "optimized model" in [6] both the proton and 
electron injection spectra are "fine-tuned" and their in- 
tensities are renormalized to explain the EGRET diffuse 
7 spectra. However, it may be not easy for the proton 
spectrum to fluctuate significantly and to be different 
from other heavy nuclei, as introduced in Q. 

It is shown that the observed peak of the diffuse 7 
spectrum at low galactic latitudes, where the dominant 
contribution is from pion decay, is at higher energies than 
the 7r° decay peak. Further the conventional model with 
reacceleration is known Q to produce less antiprotons at 
~ 2 GeV than the measurement at BESS [1] by a factor 
of about 2. Positron data also show some "excess" at 
higher energies Q. These discrepancies may all indicate 
a contribution from "exotic" sources, e.g. DMA (Toj . 
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de Boer et al. [Ill pointed out that the "GeV ex- 
cess" could be explained by the long-awaited signal of 
DMA from the Galactic halo. By fitting both the back- 
ground spectrum from cosmic nucleon collisions and the 
signal spectrum from neutralino, the lightest supersym- 
metric particle, annihilation they found the EGRET data 
could be well explained in all directions. From the spa- 
tial distribution of the diffuse 7-ray emission they con- 
structed the DM profile, with two rings supplemented on 
the smooth halo. A direction independent "boost factor" 
to the signal flux usually at the order of 100 is neces- 
sary to explain the 7-ray excess. Another factor between 
1/2 — 2 for the background flux is also needed to account 
for the spectra at different directions. However, de Boer's 
model with ring proflles and a large boost factor will lead 
to possible conflict with the antiproton flux, as shown by 
Bergstrom et al. [ll| . 

Based on the model-fitting by de Boer et al. [ll| and 
Strong's work , we try to explain the diffuse 7-ray spec- 
trum in this work by directly calculating the background 
and DMA fluxes and to overcome their shortcomings at 
the same time. By adjusting the propagation parame- 
ters we try to give consistent descriptions to the mea- 
sured spectra without any arbitrary normalization of the 
background contribution. We calculate the DMA in the 
frame of the minimal supersymmetric extension of the 
standard model (MSSM). After taking into account the 
enhancement by the existence of subhalos [U] we do not 
need the "boost factor" any more. Furthermore in our 
propagation model, we found the antiproton flux is in 
agreement with the measurements. The crucial point is 
that the enhancement by subhalos is spatial dependent in 
the Galactic halo, not "universal" as the previous works 
adopted. So the enhancement of 7-ray is different from 
that of antiproton flux, because the whole halo will con- 
tribute to the diffuse 7-ray intensity, while only antipro- 
tons produced within the diffusion region will contribute 
to the observed flux. It is found that the same scenario 
with large boost by subhalos can be used to explain the 
positron excess |14| . 
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The fluxes of DMA products are determined by two in- 
dependent factors. The first factor is related to the anni- 
hilation cross section and determined by particle physics 
of DM, while the other one is connected with the spa- 
tial distribution of DM and determined by astrophysics 
[lol |. We use the package DarkSUSY [l5| to calculate 
the particle physical factor of DMA. Scanning the pa- 
rameter space of MSSM we find the 7-ray spectrum with 
TO^ = 40 ~ 50 GeV can fit the EGRET data well. The 
branching ratios between neutralino annihilation into p 
and 7-rays are also calculated for different MSSM param- 
eters and are found to be 1/20 ~ 1/10 in a wide mass 
range. We chose a = 48.8 GeV model which predicts 
nh^ = 0.09 and g*"!^^^^! 0.055 for energies above the 
threshold Eth — 0.5 GeV. The second factor determin- 
ing the annihilation fluxes is deflned as = J jpdV 
with D the distance to the source of 7-ray production, p 
the density profile of DM and V the volume of annihi- 
lation taking place. When we consider the contribution 
from subhalos, the factor is given by the number integral 
along a direction {9,q^), <i>sub = /,^^ $"^*™d7V,„fa(0, 0). 
We use the simulation result of the subhalo distri- 
bution with mass rrigub a-t the radius r [TgI . [t^ as 
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, where 



My « 1.0 X IO^^Mq is the mass of the Galaxy, rn = 
0.14r„ ~ 36 kpc (r„ w 260 kpc is the virial radius of 
the Galaxy halo) is the core radius for the distribution 
of subhalos, r is the distance to the Galactic center (GC) 
and Nq is the normalization factor. The minimal sub- 
halos can be as light as 10~^Mq as show n by the recent 
simulation conducted by Diemand et al. [T^, while the 
maximal mass of substructures is taken to be O.OIM^ 
[1^]. The tidal effects are taken into account under the 
"tidal approximation" [16] so that the subhalos are dis- 
rupted near the GC. The total signal flux comes from the 
annihilation in the subhalos and the smooth component. 

The DM density profile within each subhalo is taken 
as the NEW [Hi, Moore [11] or a cuspier form ^ as 
P " (r/r,)-^(i+'-/r,)3-T ^ith 7 = 1.7. The last form is 
favored by the simulation conducted by Reed et al. [2l| , 
which shows that 7 = 1.4 — 0.08 log(M/M*) increases for 
smaller subhalos. We take 7 = 1.7 for the whole range 
of subhalo masses as a simple approximation. The small 
halos with large 7 1.5 ~ 2 are also found by Die- 
mand et al. [l3|. To determine the profile parameters, 
we also need to know the concentration Cy as a function 
of halo mass. Here we adopted the semi-analytic model 
of Bullock et al. [2^], which describes Cy as a function of 
virial mass and redshift. We adopt the mean Cy — rrigub 
relation at redshift zero (see also Eig. 1 of Ref. [ij). 
The scale radius is then determined as r"-'^™ — ry/cy, 
^nioore ^ r^f^ /0.63 Or r] = r^'^"'/(2 - 7). Another fac- 
tor determining the 7-ray flux is the core radius, rcorc, 
within which the DM density should be kept constant 
due to the balance between the annihilation rate and the 
infalling rate of DM particles [2^ . The core radius rcorc 
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FIG. 1: The astrophysics factor 'I>"''*''° (in unit of 
[GeV/cm^f kpc Sr"^) fr om different directions. The almost 
horizontal lines correspond to the contributions from subhalos 
only. 

is approximately in the range 10"^ 10^^ kpc for the 
7 = 1.7 proflle and 10~^ ^ 10~^ kpc for the Moore pro- 
file. In Fig. [1] we show the factor from the smooth 
component, the subhalos and the total contribution as 
a function of the direction to the GC. The (1)'"'*''° from 
subhalos is almost isotropic to different directions, this 
is because the DM distribution is almost spherical sym- 
metric and the Sun is near the GC. We can see that the 
largest enhancement for 7 = 1.7 subhalos at large angles 
can reach 2 orders of magnitude and depends on the value 
of rcorc, while for the Moore profile the enhancement is 
about one order of magnitude and for NEW profile only 
about 3 times larger. The $"''*''° for Moore and NEW 
profiles is not sensitive to rcore [3- We also notice that 
near the GG there is no enhancement. This is actually 
a very important difference from the model given by de 
Boer flTj where the "boost factor" is universal. Given 
the factor (I>''''*''° and the SUSY model we can predict 
the 7-ray fiux by neutralino annihilation. 

We now turn to the calculation of the background 
diffuse 7-ray emission, which consists of several com- 
ponents: the neutral pion decay produced by energetic 
interactions of nuclei with interstellar gas, emission by 
electrons inverse Compton scattering off the interstellar 
radiation field, the bremsstrahlung of electrons in inter- 
stellar medium, and the extragalactic background. We 
calculate the background diffuse 7-rays using the pack- 
age GALPROP '2^ which uses the realistic distributions 
for the interstellar gas and radiation fields and solves the 
diffusion equations numerically. 

We have paid extreme effort to calculate the back- 
ground so that we can give good description to the 
EGRET data after adding the DMA component. The 
injection spectra of protons and heavier nuclei are as- 
sumed to have the same power-law form in rigidity. We 
include the nuclei up to Z = 28 and relevant isotopes. For 
propagation, we use the diffusion reacceleration model 
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FIG. 2: B/C and proton spectrum in the present model. Lower curve for B/C is LIS, upper is the modulated, while the lower 
curve for proton spectrum is the modulated and upper is LIS. For the experimental data, see j^. 



[2^. The diffusion halo height of the propagation is 
taken as Zh = 1.5 kpc, which is different from 4 kpc as 
adopted in [ll| . A smaller Zf^ can effectively lower the 
p flux since it is only p from DMA in the diffusion region 
that can contribute to the flux observed on the Earth. 
The propagation parameters have been tuned to fit the 
B/C ratio and the local proton (and electron) spectra, 
as shown in Fig. O A major uncertainty in the models 
of diffuse Galactic 7-ray emission is the distribution of 
molecular hydrogen for the derivation of H2 density from 
the CO data is problematic ^26^. For example, the scal- 
ing factor Xco from COBE/DIRBE studies by Sodroski 
et al. [131 is about 2 — 5 times greater than the value 
given by IBoselli et al. [1^ in different Galactocentric ra- 
dius based on the measurement of Galactic metallicity 
gradient and the inverse dependence of Xqq on metal- 
licity, which is normalized to the 7-ray data [26[. An 
analysis of EGRET diffuse 7-ray emission yields a con- 
stant Xco = (1-9 ± 0.2) X 1020cm-V(K km s"!) for 
— 0.1 — 10 GeV [2^. Observations of particular lo- 
cal clouds yield lower values ^co = 0.9 — 1.65 x 10^° 
cm-V(K km s'^). Since the fit to the EGRET data 
for =0.1 — 10 GeV in 29] assumes only the back- 
ground contributions, we expect they give larger Xco 
than the case with new components, such as the consid- 
eration here. We find a smaller Xco = 0.6 1.0 x lO'^'^ 
molecules cm~^/(K km s^^) can give much better fit to 
the EGRET data below 1 GeV. We take Xco a constant 
independent of the radius R. As shown in Ref. [26| the 
simple form is compensated by an appropriate form of the 
CR sources. We have taken the radial distribution of CR 
sources in the form of (r/ro)"e~'''^'"~''°-'/''° with a = 1.35, 
f3 = 2.7, To = 8.5 kpc, and limiting the sources within 
rmax — 15 kpc, which are adjusted to best describe the 
diffuse 7-ray spectrum. 

The results are shown in Fig. [3] for six different sky 
regions as defined in [6!]. It should be noted that includ- 



ing the enhancement by subhalos dose not exclude the 
ring- like structures proposed by de Boer ,11]. That is 
natural since taking the subhalos into account only en- 
hances the signals coming from the smooth component 
but does not mimic the ring-like structure^which can fit 
the EGRET data at diflferent directions [ll|. Actually 
the ring-like structure, such as the tidal stream of dwarf 
galaxies are not unusual in N-body simulations. Observa- 
tions and simulations support such an idea that the ring 
at ^^14 kpc is from the tidal disruption of the Canis Ma- 
jor dwarf galaxy [13]. Recent result of the rotation curve 
also predicts a ring like structure at the similar position 
^31*]. From Fig. H we can see that the EGRET spectra 
in all regions are in good agreement with the theoretical 
values. It should also be noted that in our work we ad- 
just the propagation parameters in GALPROP and do 
not need an arbitrary normalization of the background 7 
rays as done in [11]. 

Finally we check the antiproton flux in this model. 
We first calculate the source term produced by neu- 
tralino annihilation adopting the same SUSY model as 
used for 7-ray calculation, ^p{r, E) = ■^^^^^^{yo('')^) 
where 4>{E) is the differential flux at energy E hy a, sin- 
gle annihilation and {p{r)^) = pl„,^^th + (pLb)- The 
contribution from the subhalos is given by (p('')s„b) — 
1^1°^ ^^■^b{m,r) {J p'^dV) ■ dm with Nsub{m,r) the 
number density of subhalos with mass m at radius r. We 
then calculate the propagation of p and its spectrum at 
Earth by incorporating the DMA signals in GALPROP. 
The propagation parameters are kept the same as the 
ones in background 7-ray calculation. 

In Fig. m we show the background, signal and total p 
fluxes in our model. The result is much smaller compared 
with [12] . Several ways are incorporated to decrease the p 
flux, while keeping 7-rays the same. The small z/j in our 
model helps to suppress the p flux from the smooth DM 
component. The contribution from the rings is found to 




FIG. 3: Spectra of diffuse 7-rays for different sky regions (top row, regions A, B, middle C, D, bottom E, F). The model 
components are tt" decay, inverse Compton, bremsstrahlung, EGRB and DMA (dark red curve). 
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FIG. 4; Flux of p after solar modulation. 

be greatly suppressed by slightly adjusting the ring pa- 
rameters: the inner ring is now located at R =3.5 kpc 
and the outer ring is moved from R = 14 kpc to 16 kpc. 
This is because the distance dependence of the propa- 
gation p i s steeper (exponential decrease) than of 
7-rays [32|- It can also be noted that the total p flux 



in the present model is still a bit higher than the best 
fit values of the observations at lower energies, however, 
it is consistent with data within la. The large error of 
the present data make it hard to give definite conclusion 
now. The future measurement from PAMELA 33] or 
AMS02 [s^l will finally determine if the present model is 
confirmed or disproved. 

In summary we calculate the Galactic diffuse 7-rays 
from CR secondaries and DMA. By building a new prop- 
agation model and taking into account the enhancement 
of DMA by subhalos the EGRET data may be explained 
without any "boost factor" . However, the ring-like struc- 
tures are still necessary. A lower Xco than previously 
used value is favored and the smaller halo height effi- 
ciently decrease the yield oip from DMA. The neutralino 
mass is in the range 40 — 50 GeV and very cuspy profile 
for subhalos are needed. The p flux coming from secon- 
daries (and tertiaries) and from DMA are consistent with 
present experimental bound in this propagation model by 
slightly adjusting the ring parameters. 
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